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We develop a theoretical framework for the dissipative propagation of quantized light in interacting 
optical media under conditions of electromagnetically induced transparency (EIT). The theory allows 
us to determine the peculiar spatiotemporal structure of the output of two complementary Rydberg- 
EIT-based light-processing modules: the recently demonstrated single-photon filter and the recently 
proposed single-photon subtractor, which, respectively, let through and absorb a single photon. In 
addition to being crucial for applications of these and other optical quantum devices, the theory 
opens the door to the study of exotic dissipative many-body dynamics of strongly interacting photons 
in nonlinear nonlocal media. 



Dissipation is commonly seen as a source of errors in 
quantum information and of undesired decoherence in 
strongly correlated many-body systems. However, recent 
work has shown that the very same mechanism can in- 
stead be turned into a powerful tool for quantum comput- 
ing jT], self-correcting quantum memories [2], continuous 
quantum communication and repeaters [31 H] , manybody 
entangled state generation [U [SHIO]. and dissipation- 
driven phase transitions [llj . A particular example, re- 
alized in recent experiments [T^HH] . is the propagation 
of quantized light fields in Rydberg media [T5HT5] un- 
der the conditions of electromagnetically induced trans- 
parency (EIT) [ 19_ . While Rydberg states provide strong 
long-range atom-atom interactions, EIT provides strong 
atom-light interactions with controlled dissipation. The 
resulting combination gives rise to strong and often dis- 
sipative photon-photon interactions [20^25] . which can 
be used to generate a variety of non-classical states of 
light [T^HHl [2^1155] and to implement photon-photon 
and atom-photon quantum gates [211 EH |39l HQ] . First 
wavefunction-based descriptions of two-photon propaga- 
tion in Rydbeg EIT media have revealed the emergence 
of correlated two-photon losses that could enable the de- 
terministic generation of single photons [13l|2T]. Yet, the 
fate of the remaining photon as well as the underlying 
dissipative many-body dynamics have remained unclear 
despite their essential role in the performance of future 
Rydberg-EIT-based nonlinear optical quantum devices. 

In this Letter, we address these outstanding questions 
and develop a theory for the dissipative many-body dy- 
namics of quantized light fields in a strongly interact- 
ing medium. In contrast to earlier studies [13 [H], our 
theory provides information about the many-body den- 
sity matrix of the light field, i.e. it faithfully describes 
the process of populating the n-photon states from the 
(n-l-l)-photon manifold as a photon scatters. In addition 
to opening the door to the study of photonic dissipative 
many-body physics, the theory allows one to compute 
the complex spatiotemporal structure of the generated 
non-classical light fields, whose understanding is crucial 
for applications. As two important examples that illus- 
trate this point and evince the power of our method, we 
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FIG. 1. Single-photon filter. A classical field with Rabi fre- 
quency SI resonantly coupling the excited state |e) to the Ry- 
dberg state |r) controls the propagation of the quantum field 
£. The EIT-compressed pulse length Lp is assumed to be 
smaller than the length L of the medium and the blockade 
radius Zb- 



consider the recently demonstrated single-photon filter 
[TJl [5T| and the recently proposed single-photon subtrac- 
tor i2(ii . In the limit of strong interactions, our approach 
yields exact solutions to the dissipative many-body dy- 
namics, provides an intuitive picture of the underlying 
physics, and highlights the importance of boundary ef- 
fects, i.e. of the entrance dynamics of the incoming pho- 
tons. These effects may be crucial for photon storage in 
a nonlinear quantum memory [12] . while the developed 
theoretical framework should enable the understanding of 
recent experiments [13) I14j beyond the limit of extremely 
weak input. 

The basic physics can be illustrated by considering the 
example of a single-photon filter shown in Fig. [T] In the 
absence of interactions, the probe field 8 couples with an 
effective spin wave of Rydberg atoms \r) to form a slow- 
light polariton ^19] . Whenever two polaritons are within 
the so-called blockade radius Z), of each other, the strong 
interactions between |r) atoms destroy EIT and lead 
to strong dissipation [21]. If the EIT-compressed pulse 
length Lp is smaller than the length L of the medium and 
the blockade radius z^, the first photon propagates with- 
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out losses under EIT conditions but causes scattering of 
all subsequent photons [21]. The density matrix of the 
first photon is obtained by tracing over all the subsequent 
photons. Since the first photon must already be inside 
the medium to cause scattering, the timing of the scat- 
tering events carries information about the first photon. 
Therefore, the transmitted single photon is impure. In 
the following, we develop a master-equation-type frame- 
work that allows one to determine the state of the output 
photons in this and other Rydberg-EIT problems, includ- 
ing those that do not satisfy the condition Lp < L, Zf,. 

Setup. — Let S^z), P^z), and S^z) be the slowly- 
varying operators for the creation of a photon, an exci- 
tation in state |e), and a Rydberg excitation |r), respec- 
tively, at position z. The operators satisfy the same-time 
commutation relations [£{z),£^z')] = [P{z), P"^ {z')] = 
[S{z), {z')] — 6{z — z'). The Heisenberg equations of 
motion inside the medium z G [0, L] are [13l [21] 1411 142] . 

dt£{z,t)^-dJ{z,t) + igP{z,t), (1) 
dtP{z, t) = -P{z, t) + igE{z, t) + mS{z, t), (2) 

dtS{z,t) ^ inP{z,t) ~ i [dz'V{z-z')S'' {z')S{z')S{z)(3) 



Here g is the collective atom-photon coupling constant, 
V{z) = Cq/z^ is the Rydberg- Rydberg interaction, time 
and frequencies were rescaled by 7 (the halfwidth of the 
|(7)-|e) transition), while z was rescaled by c/7. In these 
units, the blockade radius is given by zt, — (Cg/ri^)^/^ 
[2T] and c = 1. Outside the medium (z ^ [0,i]), S{z,t) 
and P{z,t) are not defined and {dt + dz)£{z,t) — 0. 

Incoming N -photon Fock state. — For simplicity, we as- 
sume that the incoming pulse is confined to a single ~ 
for simplicity, real - spatiotemporal mode h(t) satisfying 
/ dth?{t) = 1. Then an incoming A^-photon Fock state - 
before entering the medium - can be written as 
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|0>, (4) 



while its full density matrix at all times has the form 



N 



(5) 



ri=0 



where pn denotes the n-excitation part. Since the en- 
vironment possesses the information about the number 
of scattered photons, there are no correlations between 
the A*" -t- 1 terms in Eq. ([S]). The Heisenberg equa- 
tions of motion (l][3) yield master-equation- type evolu- 
tion equations for the matrix elements of p„. If all 
but at most one photon are scattered, the only nonva- 
cuum matrix element that survives in the output field is 
ee{x,y,t) — tT[pi{t)£^x)£{y)]. As shown in the supple- 
mentary material [43j , the dissipative propagation can be 
solved analytically for arbitrary photon number N under 
conditions of perfect EIT with Lp < L, Zf, and numeri- 
cally for N = 2 without any restriction on the experi- 



mental parameters. In the former case, the resulting dy- 
namics can be derived within a more general and simpler 
framework outlined below. 

Perfect EIT with Lp < L requires a large optical depth 
of the medium [19^ , implying that the absorption length 
is much smaller than the blockade radius Zf, and the com- 
pressed pulse length Lp. Since Lp < zi,, at most one pho- 
ton can propagate through the medium without losses. 
Then a fundamental question directly relevant to the ex- 
periments in Refs. [13l [M] is whether all N incoming 
photons are lost as they blockade each other's propaga- 
tion or whether one photon indeed survives. 

To answer this question, we work in the Schrodinger 
picture [44] and rewrite the input pulse [Eq. Q] outside 
the medium as 
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|0),(6) 



where the photons are now time-ordered. While the N 
incoming photons are in the same spatial mode and hence 
indistinguishable, the possibility of time ordering is the 
crucial conceptual step in the derivation. As the first pho- 
ton (z = 1) enters the medium, it turns into a Rydberg 
spin-wave excitation S moving at the EIT group velocity 
Vg = {n/g)'^c ^ c. Since Lp < Zf,, this single Rydberg 
excitation turns the entire medium seen by the remain- 
ing — 1 photons into a resonant two-level medium. As 
the absorption length is much smaller than Lp, all the re- 
maining photons get scattered into some other mode Q as 
soon as they enter the medium. We will later trace over 
those loss channels, so we can assume without loss of gen- 
erality that Q is also a one-dimensional mode with com- 
mutation relation [Q{z),Q^ {z')] — S{z — z') [45]. Once 
the entire pulse is inside the medium, we, therefore, have 



tjV>--->t2 
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where 



- -V^9 / dtMti)sHvgit - tl))|0). (8) 



(9) 



Tracing over Q, we obtain 
pit)^ I dxdyss{x,y,t)S\ym{Q\S{x) 



dt2N{N -l)h^{t2) 



n N-2 



h^{T)dT 



where ss{x, y, t) — (f>{x/vg — t, y/vg — t)/vg is the density 
matrix of the remaining spin wave with 



{x,y) = Nh{-x)h{-y) 



dzh^i-z) 



N-l 



,(11) 



which together with Eq. (|9| yields the dynamics inside 
the medium. For the purposes of this derivation, we have 



3 




FIG. 2. The narrowing and advancing of the produced single- 
photon intensity profile ee{x, x) with increasing input inten- 
sity. The input mode h^{x) oc exp[— a;'^/(2CT'^)] is shaded, (a) 
Input is a Fock state with the indicated value of A'', (b) Input 
is a coherent state with the indicated value of (n). 



ignored the small photonic component ee = Vgss. This 
solution has a simple physical interpretation: The trace 
of the integrand in Eq. ( |Io| ) is the probability that the 
second photon enters the medium (and immediately scat- 
ters) in the time interval [t2, t2 + dt2], while |'0t2(^)) is the 
unnormalized spin wave that would be propagating in the 
medium had we detected that scattering event. 

Transforming to a moving frame of reference, the den- 
sity matrix of the output photon becomes ee{x, y) — 
(j>{x,y) [see Eq. (11)]. This result shows that exactly one 



photon indeed survives the dissipative entrance dynam- 
ics: tr[p] = / dx(j){x,x) = 1. It also yields a remarkably 
simple result for the purity of the created photon: 



tr[p2 



N 



2N -1 



(12) 



As expected, the purity is smaller than unity because the 
timing of the scattering event carries some information 
about the remaining spin wave |V't2(i))- Crucially for 
applications, the purity does not vanish but approaches 
1/2 as N — oo. Surprisingly, it is independent of the 
mode shape h{t). Furthermore, the eigenvalues pt and 
eigenvectors (f)i{x) of (j){x,y) can be easily found j 43l by 
using the change of variables x — > dzh^ {— z) , which 
makes the density matrix and hence Pi independent of 
h{t). Physically, this surprising behavior emphasizes the 
fact that the key role is played simply by the arrival order 
of N identical photons and not by the shape of the mode. 

This dynamics at the medium boundary leads to a 
slight narrowing and advancing of the single-photon pulse 
(f>{x,x), as shown in Fig. [2][a) for a typical Gaussian in- 
put mode. This behavior can be traced back to the first 
scattering event, which projects the leading photon into 
the medium. This effect becomes more pronounced with 
increasing N, since the larger N is the sooner the first 
scattering event takes place. More succinctly, the prob- 
ability distribution of the first photon is obviously ad- 
vanced and narrower relative to the normalized probabil- 
ity distribution h'^{t) of the entire incident pulse. Fortu- 
nately, 4>(x,x) and (jui^) shorten extremely slowly with 
iV as ^ 1/ \^\og N, keeping the associated losses at a min- 
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FIG. 3. (a) For an incoming N — 2 Fock state, the two- 
photon amplitude and the single-spin-wave density matrix at 
different times, (b) Approximate analytical (dashed) and ex- 
act numerical (solid) results for the efficiency tr[pi] (red) and 
purity tr[pi]/tr[pi]^ (blue) of the single excitation. The pur- 
ple arrow indicates the efficiency of retrieving Eq. ([9| with 
N = 2. (c) As a function of (n) of the incoming coherent 
pulse, for the indicated values of the blockaded optical depth 
ODb, the estimated efficiency rj of the single-photon source. 



To verify this intuitive picture, we carried out nu- 
merical simulations for N = 2 incoming photons and 
h{t) cx 1 - 4(i/r - 0.5)2 g [0, T]) using the full propa- 
gation equations derived from Eqs. (l][3) [13] ■ This form 
of h{t) is motivated by optimal photon storage [i^ I46j . 
The results are shown in Fig.[3]ja,b). In Fig.[3](a), the left- 
bottom quadrant corresponds to both photons being still 
outside the medium, so t = is described by Eqs. (4|6|. 
The top-right quadrant corresponds to both photons be- 
ing inside the medium, so t = T is described by Eq. (|9|. 
Finally the remaining two quadrants correspond to the 
first photon being already inside the medium while the 
second photon is still outside. Fig. [3f^b) shows a compar- 
ison to the analytical prediction from Eqs. ( 9|11 ) |47j . 
While imperfections keep the single-photon conversion 
efficiency slightly away from unity, the overall physical 
picture is very well confirmed by our numerical simula- 
tions. To verify that losses induced by the finite width of 
the EIT transparency window - and not the correlated 
photon dynamics upon pulse entrance - constitute the 
dominant imperfection, the purple arrow in Fig. [3](b) in- 
dicates the efficiency of retrieving Eq. ^ with N = 2 
back out of the medium. 

Arbitrary incoming state. — Since any mixed state can 
be represented as a classical mixture of pure states, it is 
sufficient to consider an arbitrary pure input state 



|^)=^c™|m). (13) 

m 

Following the above treatment, we trace over all photons 
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except for the first one to obtain 

p - (co|0) + ci|l))(cS(0| + cl(l|) + |c„|2p„,(l4) 

m>2 

where pm is the single photon obtained from \m). The 
single-photon conversion efficiency is 1 — |cop, i.e. lim- 
ited by the vacuum component of the input state. The 
corresponding purity is tr[p^]/(l — |coP)^, where 

tr[p?]= E I^-H^-I't^^T^^^^^T^^- (15) 

m,n> 1 



(n- 



l)(n- 



For a coherent input, |cip = e ^"^(rt)Yz!, the efh- 
ciency is thus 1 — e^^"\ while the single-photon purity 
is (1 - e-<">)-2(i - e-2(">(i + 2(n)))/2, which falls off 
monotonically from 1 to 1/2 with increasing (n). Since 
I Cop drops exponentially with (n), a small average num- 
ber of incoming photons (n) 10 is sufficient to make the 
single-photon source deterministic. Repeating the above 
derivations, one obtains for the output density matrix 



{x,y) = {n)h{~x)h{-y)exp 



(n) h'^{-z)dz 

J min(a:,y) 



(16) 



which can easily be diagonalized [13]. As for a Fock- 
state input with Gaussian h{x), the output pulse shortens 
extremely slowly with increasing (n) (~ 1/ yJ\og{n)) , as 
shown in Fig. [2]jb) . 

The efficiency of this single-photon source - imperfect 
due to the finite width of the EIT window - can be es- 
timated from the analytical form of the density matrix 
without involving interactions. We assume that the in- 
coming pulse is stored without interactions into the spin- 
wave s{z) C!C 1 — 4:{z/L — 0.5)^ and that the single-spin- 



wave density matrix [Eq. (11) or Eq. (16l] is retrieved 
forward. The efficiency rj of the single-photon source can 
then be estimated as the product of these two - storage 
and retrieval - efficiencies. The (n)-dependence of r/ for 
a coherent input is shown in Fig. ptc) for different block- 
aded optical depths ODh, assuming the entire medium is 
blockaded. The relatively poor scaling of the efficiency 
with ODf, results from the cusp of the density matrix 
4>{x,y) along the diagonal {x = y), which carries high- 
frequency components. In a magneto-optical trap (den- 
sity Af ~ 10^2 cm-3), ODb - 10 [13], and hence r] w 0.2 
can be achieved. In a BEG [371 iS], AT ~ lO" cm^^ ^^^-^ 
give ODb ~ 1000 and rj « 0.9. The efficiency can be 
further increased by using photonic waveguides [i5ti52] 
and by further optimizing h{t) and retrieving backwards 

[Hill. 

Despite their impurity, the single photons produced 
with this method are a valuable resource. In particular, 
the impurity would not interfere with applications that 
do not rely on Hong-Ou-Mandel interference [S3], such 
as optical quantum computing with impurity-insensitive 
two-qubit gates (e.g. [54]) or the BB84 quantum key dis- 
tribution protocol |55j . For applications that rely on 
Hong-Ou-Mandel interference, the photon can be puri- 
fied in the following ways. First, the detection of the 



first scattered photon would yield a pure photon. Sec- 
ond, the impure single photon can be purified to the 
dominant eigenvector with probability pi (pi = 0.69 
for (n) ^ 1). This can be accomplished, e.g., using an 
atomic ensemble in a cavity to store only the mode 
\(pi), heralded by the absence of a click at the cavity out- 
put, followed by retrieval, which can be shaped into any 
desired mode |57j . 

Photon subtraction. — To demonstrate the versatility of 
the developed theory, we now apply it to the Rydberg- 
based single-photon subtractor proposed in Ref. |26j . 
showing that this scheme also yields impure output 
states. The detailed physics of such a setting [43J is com- 
plementary to the single-photon source in so far as the 
density matrix of the remaining photons is obtained by 
tracing out the first one. Since the timing of the absorp- 
tion carries information about the remaining photons, 
the density matrix of the latter is impure as well. In 
fact, the single-photon subtractor and the single-photon 
filter complement each other to make the original pure 
state. Hence, the impurity and the entire eigenspectrum 
of the reduced density matrix are identical in the two 
cases. 

This can be shown by tracing over the first photon in 
Eq. ( 13 ) to obtain 



where 

m>l 



dt,h^ti)\iJtAt)){i^tAt)\, (17) 



dt'h{t')£\t-t') 



|0>, 



which has the same eigenspectrum as Eqs. (14 1. 

Outlook. — In conclusion, we extended the dynamics of 
open quantum systems of Rydberg atoms [551 - 157] to in- 



clude the dissipative quantum dynamics of the propa- 
gating light field, which is crucial for the understand- 
ing of recent experiments [T31 [H] . While we focused on 
the case of blockade radius that is larger than the size 
of the EIT-compressed pulse, the developed framework 
also applies to finite blockade radii [lU [14] and can be 
easily extended to a time-dependent blockade radius, as 
relevant for photon storage via time-dependent control 
fields. Extensions to non-dissipative unitary evolution 
|40j . media with longitudinal density variations, incom- 
plete transverse blockade, as well as finite Rydberg-state 
lifetime are straightforward 13 . Finally, we expect our 
calculations to be extendable to other light-processing 
modules, such as the quantum filter of Refs. [68j 69]. 
Most importantly, our approach may lead to a simplified 
effective theory for the many-body dissipative dynamics 
of correlated photons in strongly interacting media. 
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I. EQUATIONS OF MOTION FOR TWO 
INCOMING PHOTONS 



In this Section, we give the equations of motion that 
were used to do the numerics for two incoming photons 
and to obtain Fig. 3(a,b) in the main text. 

In the case of two incoming photons, the full density 
matrix 

Pit) = e(t)|o)(o| + + \Mt)){Mt)\ (SI) 

consists of the unnormalized two-excitation wavefunction 
\Mt)) = \hx I dyEE{x,y,t)£\x)£\ym 



+ I dxj dyEP{x,y,t)E\x)P\y)\Q) 
dx I dyES{x,y,t)8\x)S\y)\Q) 



+ i / dx j dyPP{x,y,t)P\x)P\vm 



dx j dyPSix,y,t)P\x)S\y)\0) 

dx J dySS{x,y,t)SHx)S\y)\0),iS2) 
the unnormalized single-excitation density matrix 
pi(t)= fdxj dyee{x,y,t)i\y)\0){0\£ix) 

dx J dyep{x,y,t)P^iy)\0){0\£{x) 

+ I dxj dype{x,y,t)£\y)\0){0\Pix) 

dx J dyesix,y,t)S\y)\0){0\i{x) 

' dx J dyseix,y,t)£\y)\0){0\S{x) 

' dx J dypp{x,y,t)PHy)\0){0\P{x) 

" dx j dyps{x,y,t)SHymmP{x) 

' dx j dysp{x,y,t)P\y)\Q){Q\S{x) 

' dx I dyss{x,y,t)S\ym{Q\S{x),{Si) 



and the vacuum component e(t)|0)(0|. Here / integrates 
over (— cx),oo), while /' integrates over [0,L]. Without 
loss of generality, we take EE, PP, and 55* to be sym- 
metric [e.g. EE{x,y) = EE{y,x)]. If the input state 
had correlations between different Fock states, one would 
need to include coherences between manifolds of different 
photon number; the method we discuss can be naturally 
generalized to these situations. 

All terms in p{t = 0) vanish except for EE{x,y,0) = 
\/2h{—x)h{—y), where we assume h{t < 0) = 0. The 
equations of motion for EE, EP, ES, PP, PS, and SS 
can be obtained by expressing them in terms of {^2} [e.g. 
ES{x,y) = {Q\£{x)S{y)\il)2{t))] and using Eqs. (1-3) in 
the main text. For a; ^ [0, L], y G [0, L], they are 



{dt + d^, + dy)EE = igEP, 
{dt + d^ + l)EP = igEE - 
{dt + dx)ES = inEP, 



iQES, 



(S4) 
(S5) 
(S6) 



and describe the EIT propagation of photon y, while pho- 
ton X propagates outside the medium with the speed 
of light (c = 1 in our units). Using EE{x,y,t) — 
\f2h(t — x)h{t — y) to set the boundary conditions at 
y — 0, these equations are solved for a; < 0, y G [0, L] 
to give the boundary conditions for the equations in the 
region x, y G [0, L]: 



{dt +d, + dy)EE 


- igEP+, 






(S7) 


{dt -t- a, + l)EP 


= i9{EE 


^PP)+i 


QES, 


(S8) 


{dt + d^)ES 


= igPS + 


iQEP, 




(S9) 


{dt + 2)PP 


= i9EP+ 


+ inES+ 




(SIO) 


{dt + l)PS 


= igES + 


in{pp + 


SS), 


(Sll) 


{dt + iV{r))SS 


= inps+ 






(S12) 



where EP±{x,y) = EP{x,y) ± EP{y,x), ES±{x,y) = 
ES{x, y) ± ES{y, x), PS±{x, y) = PS{x, y) ± PS{y, x), 
and r — X ~ y. The solution to these equations can then 
be used to set the boundary conditions at a; = L for Eqs. 
( p4p6| in the region x > L, y Cz [0, L], which can in turn 



be used calculate the outgoing two-photon pulse. 

Now we turn to the evolution equations for the single- 
excitation density matrix pi. We first note that 



es{x,y) = {£^{x)S{y)) - J dzEE* {x, z)ES{z,y) (S13) 
-j dzEP*{x,z)PS{z,y)- j dzES* {x, z)SS{z,y). 
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The equation of motion for {E^^ (x) S (y)) follows from Eqs. 
(1-3) in the main text. Together with the equations of 
motion for the two-photon amplitudes, this yields the 
equation of motion for es(x,y), and, similarly, for all 
matrix elements of pi. The following source terms will 
describe the transfer of population from \ip2) to pi: 



fee{x,y) ^2 J 
fep{x,y) = 2 
fes{x,y) = 2 



fpp{x,y) = 2 



fpsix,y) 



dzEP*{x,z)EP{y,z), (S14) 

dzEP*{x,z)PP{y,z), (S15) 

dzEP*{x,z)PS{z,y), (S16) 

dzPP*{x,z)PP{y,z), (S17) 

dzPP*{x,z)PS{z,y), (S18) 



fss{x,y)^2 / dzPS*{z,x)PS{z,y). 



(S19) 



As expected, in the interaction-free case {V — 0) and 
assuming perfect EIT, the source terms vanish because 
|e) is never populated, so all components of IV'2) involving 
P vanish. With these definitions, for x,y G [0, L], 

[dt + dj;+ dy)ee = ig{ep - pe) + fee, (S20) 
(dt +dj: + l)ep = ig{ee - pp) + iUes + fep-, (S21) 
{dt + dx)es = iQep-igps + fes, (S22) 
{dt + 2)pp = ig{pe - ep) + iVL{ps - sp) + fpplS23) 
{dt + l)ps = -iges + in{pp - ss) + /p^, (S24) 
dtss ^ iVt{sp- ps) + fss, (S25) 

while pe{x,y) = ep*{y,x), se{x,y) = es*{y,x), and 
sp{x,y) = ps*{y,x). Equations of motion outside of 
x,y € [0, L] can be obtained in the same way. 

We do numerical calculations in the regime of good 
EIT (< 1% single-photon loss). Thus, to a good approx- 
imation, photon scattering occurs only when both pho- 



tons are inside the medium. We therefore solve Eqs. ( S20 



S25 1 with vanishing initial and boundary conditions. 

We note that the equations can easily be extended [SI' 
to include longitudinally varying density, finite decoher- 
ence rate of S, as well as cases where the blockade radius 
is smaller than the transverse extent of the probe beam. 



II. IDEAL SINGLE-PHOTON GENERATION 
FROM 2 PHOTONS 

In this Section, in the case of an input Fock state with 
= 2, we show how Eqs. (9,11) in the main text arise 

from the full equations of motion presented above. 
Let's assume that EIT is perfect and that Lp < L < Zf,. 

Then, for a; < and y e [0, L], Eqs. ([S4][S6]) give 



where Vg = {ft/g)'^ in our units. From Eqs. ( S9|S11 ), we 
obtain dtES « -d^ES - g'^ES, so that, for x~ye~[0, L], 



0,y,t 



PS{x, y, t) w igES{x, y, t) igES{x 

« ~ig.f2/7gh{t)h{t - y/vg)e'a 



X e 



-g X 



(S27) 



which describes the absorption of the two-excitation am- 
plitude over the absorpti on le ngth 1/g^ <C Lp. Inserting 
this expression into Eq. ( S19[ ), we obtain 



fss{x, y, t) = 2h\t)h{t - x/vg)h{t - y/vg)/vg.{S28) 



Then from Eqs. ( |S22|S24|S25 1, to a good approximation, 

(S29) 
(S30) 



es = —{V,/g)ss + {Q/ g^)dxSS, 



dfss = ~2n^ 



Qg{se + es) + f^. 



Inserting Eq. (S29l into Eq. (S30), we obtain 



dfSs = -VgdRss + fss{x,y,t), 



(S31) 



which describes how the source fss puts excitations into 
ss; and, as soon as the excitations are put in, they start 
moving at Vg along R = {x+y)/2. This equation assumes 
the wavepacket's frequency components all fit inside the 
EIT transparency window. Solving this equation gives 



ss{x,y,t) 



-h{t-—)h{t~^) 



dt'h\t')[S32) 



which, for N — 2, generalizes Eqs. (9,11) in the main 
text to cases when the pulse has only partially entered 
the medium. Eq. ([S32| yields Tr[p2]/Tr[pi]2 = 2/3 for 

all t and Tr[pi] = dTh'^{T) , which give the dashed 

lines in Fig. 3(b) in the main text. To a good approx- 
imation, pi satisfies the dark-state-polariton condition 
= ~\f^^^ ~ 'VgSs. This derivation can easily be ex- 
tended to include the effect of a finite EIT transparency 
window width, which partially explains the slight dis- 
crepancy between analytics and numerics in Fig. 3(b) in 
the main text. 



HI. IDEAL SINGLE-PHOTON GENERATION 
FROM iV PHOTONS 

In this Section, we generalize Sec. |TT]to arbitrary N . 

Let Xm = xi, . . . , Xra, Era = E . . . E (where E is re- 
peated m times to denote the m-photon wavefunction) , 
and h{t — x,„) = YViLi ^(^ ~ ^i)- Then, for x^r < 0, the 
incoming A^-photon state is given by 



EN{^N) = VN^.h{t-^N). 



(S33) 



Once the first two photons enter the medium (xAr_2 < 



ES{x, y, t) = ^^2/^gh{t - x)h{t ~ y/vg), (S26) = -ig^N\/vgh{t - XN-2)h{t)h{t - XN/vg)e-^' 



and xn-i,xn > 0), we have, in analogy with Eq. (S27l, 
En-2PS{^n) (S34) 
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So, by analogy with Eqs. ( |S19|S28| , 

/ejv-2seAr-2s(^-'V-li ^W-l) (S35) 
= 2jdzEN-2PS* {xN-2, Z,XN-i)En-2PS{-X.'n_2, Z,x'ff_^) 



h(t-XN_2)h(t-x.L_2)h^(t)h(t 

Vn V. 



Applying group ve locity propagation along (^jv-i 
x'^_i)/2 [as in Eq. (fSSll)], we have [as in Eq. (IS35|] 



iV! 

ejv-2se7v-2s(xAr_i, x^_i) = — h{t - XN-2)h{t - xjv-2) 

Vn 



Vg Vg Jf. ""'H^W-l>^iV_l) 

^'9 



Allowing now the third photon to enter the medium 
{xn-2,x'^-2 > 0), we have 

(X7V-2,X^_2) (S37) 
= eAr_2SeAr_2s(xAr_3, 0, a;jv-2, x'^_3, 0, 

m. 



-h{t - XN-3)h(t ~ x'j^_3)h\t) 



XN-2 



dt'h^t'). 



Applying group velocity propagation along {xn^2 H 
x'^_2)/2, we have 

eAr_3SeAr„3s(xAr_2, x'^_2) 

= —h{t - xjv-3)^(i - x'^r_3)/^(^ - -^^)h{t 

Vn Va Vn 



min(x„_2,x^_2) 



dt'h^{t') 



(S38) 



Allowing the fourth photon to enter the medium, we have 

/ejv-4seN-4s(xAf-3,X^„3) (S39) 
= eAr_3SeAr„3s(xAr_4, 0, a;Ar_3, x'^„4, 0, a;^_3) 
iV!. 



'-h{t - y.N-i)h{t ~ y.'pf_^)h [t) 



xh{t )h{t )- 

Vn Vn 2 



Applying group velocity propagation along (a;Ar_3 + 
x'^_^)/2, we have 



eAr_4SeAr_4s(xAr_3, x'^_3) 

= — '-h{t — xi) . . . h{t — Xjq_^)h[t ^^—^)h(t 



Vn 



3! 



^ dt'h'it') 



min(xjv_3,xjy_3 



(S40) 



We continue in this way until we reach 

ss(xi,Xi) = 

= —h{t -)h{t ^) 



dt'h'^it') 



-^ N-1 



which generalizes Eqs. (9,11) in the main text to cases 
when the pulse has only partially entered the medium. 



IV. EIGENVECTORS OF THE 
SINGLE-PHOTON DENSITY MATRIX 



In this Section, we study the eigenvectors and eigen- 
values of the single-photon density matrix, Eqs. (11) and 
(16) in the main text, obtained via single-photon filtering 
from Fock-state and coherent-state inputs, respectively. 

We first study the eigenvectors and eigenvalues Pi of 
the single-photon density matrix (j){x, y) given in Eq. (11) 



dzh^ 



we ob- 



in the main text. Defining x = J_ 
tain p — Jq dxdy4>{x,y)£\y)\0){0\£{x), where 4'ix,y) = 

N[mm{i,y)f-\ £{i) = £{x)/h{-x), [£{x),t{y)] = 
5{x — y). The eigenvalues Pi are then the solutions of 

the characteristic equation J^i/n 2-\/ {N — 1)/ {Np) — 

0. In particular, in the limit N ^ oo, pi are the 
roots of Jo[2/y/p] — 0. The eigenvectors of (l){x,y) are 

Mi) oc Ji_i/Ar [2^iN - l)/{Np,)i^/^] . In 

particular, for N — 2, pi — 27r^^ (n — ^, 0i(i) = 
-\/2 sin [vr (n — |) 5] ■ While (j){x, x) and 0i(i) shorten as 

1/N with increasing TV, 4>{x,x) and 0i(a;) — h{—x)4>i (i) 
shorten much slower as 1 / \/\og N for a Gaussian h{x). 

We now study the eigenvectors of the single-photon 
density matrix (f)(x,y) given in Eq. (16) in the main 
text. Following the same change of variables, we obtain 
(f>{x,y) — (n) exp [— (n) (1 — min(i;, y))], which, for {n) ^ 

1, agrees with the above (j){x, y) provided one identifies N 
with (n). For general (n), the eigenstates (pi of 0(i;, y) are 
linear combinations of e-<"Hi-i)/2 ji[2e-<"Hi-5)/2/^] 

and e-<">(l-*)/2Yl[2e-<">(l-^)/V^/P^], where p^ are the 
eigenvalues. 



V. SINGLE-PHOTON SUBTRACTION 

In this Section, we present a formal derivation of Eq. 
(17) in the main text, which describes the output of a 
single-photon subtractor |S2] . In addition to verifying 
Eq. (17), this method allows one to treat deviations from 
the ideal result. 

Following Ref. |S2j . the atoms can be in one of two 
collective states \G) and \E). The density matrix then 
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evolves according to the following master equation: 

1*00 

f>^-i{Ho,p]+V / dx 2£{x)\E){G\p\G){Ep{x) 
Jo L 

-£^x)S{x)\G){G\p - p\G){G\£\x)S{x)\ . (S42) 



Hq here describes simple propagation of hght in vacuum. 
The photon is subtracted within a few absorption lengths 

of x = 0, so the remainder of the medium plays no 
role provided Zb > L; hence we assumed L — 00. 

Here, for simpHcity, we only present the derivation for 
two incoming photons |2). GeneraUzation to an arbitrary 
incoming state |^) = J2n ^^n\n) is straightforward. 

Therefore, the full density matrix 

P = Pi + \^2){H (S43) 
consists of the two-photon wavefunction 

IV2) = IJ dxdyEE{x,y)£\x)£\ym\G) (S44) 

and of the single-photon density matrix 

pi= j dxdyee{x,y)£^{y)mE){EmS{x). (S45) 



One then finds the following equations of motion: 
dtEE{x,y) = -d^EE - dyEE - T[H{x) + H{y)]EE, 

dtee{x,y) = -d^ee - dyce + 2T dzEE{y,z)EE*{x,z), 

Jq 

where H{x) is the Heaviside step function. 

Starting with the boundary conditions EE{x,y,t) = 
\/2h{t — x)h{t — y) for x,y <0, we solve for EE: 

EE{x, y, t) = s/2h{t - x)h{t - y)e-^[^('^)^+^(^)^l(S46) 

Inserting this into the equation of motion for ee and using 
the fact that the absorption length is much shorter than 
the (now uncompressed) pulse duration, we obtain 

dtee(x, y, t) = —dxee — dyce 

+2h\t)hit - x)h{t -y)[l- H{x)][l - H{y)].iS47) 
For x,y < 0, this can be integrated to give 

ee{x,y,t) = 2h{t-x)h{t-y) f h^{t'), (S48) 

J — 00 

SO that, in the remaining three quadrants of the xy plane, 

/t— in£ix(a;,y) 
rft'/i2(i')(S49) 
-00 

which is a special case of Eq. (17) in the main text. 
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